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ABSTRACT 

We consider the present-day photometric and chemical properties of elliptical galaxies, 
adopting the bimodal star formation scenario of Elbaz, Arnaud & Vangioni-Flam (1995). 
These models utilise an initial mass function (IMF) biased heavily toward massive stars 
during the early phases of galactic evolution, leading to early Type II supernovae- 
driven galactic winds. A subsequent lengthy, milder, star formation phase with a normal 
IMF ensues, responsible supposedly for the stellar population observed today. Based 
upon chemical evolution arguments alone, this scenario has been invoked to explain the 
observed metal mass, and their abundance ratios, in the intracluster medium of galaxy 
clusters. Building upon the recent compilations of metallicity-dependent isochrones for 
simple stellar populations, we have constructed a coupled photometric and chemical 
evolution package for composite stellar populations in order to quantify the effects of 
such a model upon the photo-chemical properties of the resultant elliptical galaxies. We 
demonstrate that these predicted properties are incompatible with those observed at 
the current epoch. 

Key words: galaxies: elliptical - galaxies: evolution - galaxies: abundances - galaxies: 
intergalactic medium - galaxies: clusters 
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^ X-ray observations of the hot, optically thin, gas compris- 
^ ing the intergalactic medium of clusters of galaxies have es- 
tablished firmly the presence of a metal enriched compo- 
nent of non-primordial origin (e.g. Sarazin 1986 and Mat- 
teucci & Vettolani 1988, and references therein). Arnaud 
et al. (1992) have shown that a reasonably tight cor- 
relation exists between the total V-band luminosity of a 
cluster's elliptical-plus-lenticular population and the clus- 
ter's intracluster medium (ICM) iron mass - in particular, 
Mpg^ Si 0.02iy'''^'' . Determining a-element to iron ratios 
for cluster ICMs is somewhat more problematic as the gas 
temperatures involved are usually too high to allow observa- 
tion of emission lines besides iron. Recent ASCA measure- 
ments for the few clusters which have fairly secure oxygen 
and silicon detections do seem to indicate that an overabun- 
dance relative to iron exists, with [O/Fe] and [Si/Fe] rang- 
ing from -1-0.20 -|-0.70 and -|-0.10 -|-0.60, respectively 
(Mushotzky 1994). 

Competing theories as to the origin of this metal en- 
riched component of the ICM tend to fall into two camps 
- (i) ram pressure stripping of galactic interstellar media 
(ISMs) (Sarazin 1979), and (ii) supernovae (SNe)-driven 



winds subsequent to the initial phase of intense star forma- 
tion in cluster ellipticals (Larson & Dinerstein 1975). Based 
upon cluster energetics and the lack of significant abun- 
dance gradients in the cores of the richest clusters. White 
(1991) has argued that the "stripping" scenario is not a vi- 
able one for most clusters. Additional evidence in support 
of the favoured "wind origin" for ICM metals can be seen 
in a number of recent papers, including Matteucci & Vet- 
tolani (1988), David, Forman & Jones (1991), Ciotti et al. 
(1991), Arnaud et al. (1992), Renzini et al. (1993), Mihara 
& Takahara (1994), Elbaz, Arnaud & Vangioni-Flam (1995), 
and Matteucci & Gibson (1995). 

We note here, and simply refer to the aforementioned 
papers for further proof, that most of these wind models 
tend to underproduce iron relative to the observed relation 
of Arnaud et al. (1992) by factors of ~ 3 ^ 30, for models 
which honour both the Salpeter (1955) solar neighbourhood 
initial mass function (IMF) and the [a/Fe]si -|-0.4 seen in 
the ICM. 

A typical "fix" for this situation has been to resort 
to an IMF for ellipticals which is biased toward massive 
stars (e.g. David, Forman & Jones 1991; Mihara & Taka- 
hara 1994; Matteucci & Gibson 1995), increasing the yield 
of Type II SNe-originating iron by suitably (and more im- 
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portantly, freely) adjusting the IMF slope, star formation 
rate, and/or galactic wind epoch (i.e. the time at which the 
bulk of star formation is taken to cease) to match the ICM 
iron mass and a-element ratios. Corroborating evidence for 
a flatter IMF in ellipticals usually references the population 
synthesis work of Arimoto & Yoshii (1987). As an aside, 
certain models either choose to relax, or were unaware of, 
the [O/Fe]>0 constraint (e.g. Matteucci & Vettolani 1988; 
Mihara & Takahara 1994) and will not be discussed further 
here. 

An alternate fix to the iron underproduction "problem" 
retains the preferred Salpeter (1955) IMF for the bulk of the 
star formation phases, but postulates the existence of an 
earlier, more intense, stage of star formation (i.e. 10^ 10^ 
Mg/yr) with a Salpeter (1955) IMF truncated at a high 
lower mass limit of 3 Mq, reminiscent to that put forth 
for starburst galaxies (e.g. Chariot et al. (1993), and ref- 
erences therein). The initial star formation terminates after 
~ 20 Myr due to the onset of a powerful Type 11 SNe-driven 
wind which expels, over the subsequent ~40 Myr, ~l/2 the 
initial galactic mass to the cluster ICM in the form of highly 
enriched gas. The second phase of star formation, with the 
"normal" IMF, then proceeds, the stars being formed out 
of the enriched material still present at the end of the wind 
phase. Again, by fine-tuning the star formation parameters 
and wind onset and duration, it is fairly straightforward to 
eject both the necessary amount of iron and in the propor- 
tions dictated by the ASCA observations. As pertains to 
elliptical galaxies and the ICM, this bimodal star formation 
model is introduced in the paper by Arnaud et al. (1992), 
and discussed in greater quantitative detail by Elbaz, Ar- 
naud & Vangioni-Flam (1995). It is to the predictions of 
this latter work to which we will concentrate in the analy- 
sis which follows, although the conclusions are valid for the 
simple bimodal model as a whole. 

The great danger in adopting either of the above "fixes" 
is that it treats cavalierly, or as is more usually the case, 
ignores entirely, the impact upon the resultant model ellip- 
ticals' photometric properties. It is of questionable use to 
construct synthetic elliptical models which explain the ICM 
abundances, but bear little resemblance to actual observed 
galaxies! 

Unfortunately, until the very recent past it was ex- 
tremely difficult to couple either of the above "chemical 
evolution motivated" models to photometric predictions due 
to the lack of non-solar metallicity stellar evolution tracks 
and photometric calibrations. It has only been with the re- 
lease of Kurucz's (1993) model stellar atmospheres that a 
reasonably complete, and accurate, grid has become avail- 
able for such an undertaking. The subsequent compilations 
of metallicity-dependent isochrones based upon the Kurucz 
grid, by Bertelli et al. (1994) and Worthey (1994), now 
make it possible to re-examine the photometric predictions 
for some of the galactic wind models adopted to explain the 
ICM abundances. To this end we have constructed a pop- 
ulation synthesis package which builds upon the isochrones 
listed above and which is coupled to our chemical evolution 
code described elsewhere (Gibson 1995). 

The bimodal star formation/IMF models mentioned 
earlier are successful at reproducing many of the chemi- 
cal properties of elliptical galaxies. Because of their extreme 
iron production during the early phase of star formation 



(and its subsequent ejection to the ICM) in their models 
(e.g. ten to twenty times greater than the Case A models 
of Mihara & Takahara (1994)), and the above "favourable" 
chemical predictions, we felt it imperative to ascertain if the 
model was equally successful in replicating the present-day 
photo-chemical properties of the underlying ellipticals. Self- 
consistent modeling of the photo-chemical properties of bi- 
modal models has never been demonstrated before. We stress 
that our goal is not to claim to have solved the so-called 
iron underproduction problem (although, see Gibson & Mat- 
teucci (1995) for the latest approach using the more con- 
ventional single IMF approach), but only to draw atten- 
tion to shortcomings in their proposed bimodal star forma- 
tion/galactic wind solution. 

Section 2 introduces the photometric evolution tech- 
nique adopted, and highlights some subtle differences be- 
tween the Worthey (1994) and Bertelli et al. (1994) metal- 
rich isochrones. Missing low-mass extensions to the latter set 
are also described. The relevant equations in the Elbaz, Ar- 
naud & Vangioni-Flam (1995) bimodal star formation sce- 
nario are then introduced, and the colours and metallici- 
ties of their favoured models presented in Section 3. Their 
predicted colour-luminosity-metallicity relation is then com- 
pared with observations, and the results summarised. 



2 ANALYSIS 

2.1 Isochrone synthesis 

As the basis for our photometric properties analysis, we 
build upon the metallicity-sensitive isochrones of Worthey 
(1994) and Bertelli et al. (1994). General photometric evo- 
lution techniques and history are discussed at length in these 
two papers, and the numerous references therein. Both com- 
pilations span ~2.5 dex in [Z] ([Z]=log Z - log Zq) and sam- 
ple the primary stellar evolutionary stages from the zero 
age main sequence to well down the white dwarf cooling 
curve (or carbon ignition, depending upon the initial mass). 
Worthey (1994) does not include any post-asymptotic giant 
branch contribution in his tables, but comparisons done us- 
ing Bertelli et al. (1994), with and without the planetary 
nebulae/white dwarf stages, showed this to be an unimpor- 
tant omission for the analysis which follows. 

As the Bertelli et al. (1994) isochrones have a lower 
limit of 0.6 Mq, it was necessary to extend these down 
to ~0.1 Mq, in order to include the contribution from M- 
dwarfs, and fully sample the full range of masses in the IMF. 
To this end, we transformed the newly available low temper- 
ature stellar models of AUard & Hauschildt (1995), and in 
particular, their log g = -\-5.0 grid, from the theoretical plane 
in which they are provided, to the observer's plane, via the 
mass-luminosity-temperature relationships shown in Figure 
1. These were derived from the older models of VandenBerg 
et al. (1983), so while not ideally matched to the AUard & 
Hauschildt (1995) models, they were adequate for our pur- 
poses. Interpolation and extrapolation for the Bertelli et al. 
(1994) isochrone metallicities was done where necessary, and 
continuity at m Si 0.6 Mq in the AfBoi-Teff and My-colour 
planes between the two sources was ensured. As in Bruzual 
& Chariot (1993), the colour properties of these low tem- 
perature extensions were assumed to be unevolving. 
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Figure 1. Effective temperature-mass-iuminosity reiations used to transform Aiiard &z Hauscliiidt (1995) iow temperature modeis to tlie 
observer's piane. 



Worthey (1994), on the other hand, used the M-giant 
fluxes of Bessell et al. (1991), supplemented with optical 
spectrophotometry from the Gunn & Stryker (1983) stel- 
lar atlas, to extend his isochrones down to m Si 0.1 Mq. 
Worthey's (1994) low temperature extension is virtually iden- 
tical to that adopted above for the Bertelli et al. (1994) 
isochrones. 

Table 1 lists the integrated V-K colours for simple stel- 
lar populations (SSPs) using the various isochrones of both 
Bertelli et al. (1994) and Worthey (1994). We have consid- 
ered other colours in the analysis which follows, but we will 
stress V-K, as it provides a more valuable constraint than, 
say, B-V. Over the range of luminosities considered in our 
study, B-V does not vary by more than ~ 0.2 magnitudes, 
whereas the V-K versus My relation is considerably flatter 
and spans > 1 magnitude. As can be seen in Figure 7b of 
Arimoto & Yoshii (1987), the near-vertical distribution of 
ellipticals in the B-V versus My plane makes this colour, by 
itself, a poorer constraint. 

We do not wish to enter into a detailed analysis of the 
differences between the two sets of isochrones. Work of this 
nature can be found in Worthey's (1994) paper, and the im- 
portant study by Chariot, Worthey & Bressan (1995), inter- 
comparing the models of Worthey (1994), Bressan, Chiosi 
& Fagotto (1994), and Bruzual & Chariot (1993). 

The one point which we do want to highlight in Table 
1 is the behaviour of the super-solar colours. It is apparent 
that for metallicities Z<Zq the colours agree to within a 
few hundredths of a magnitude. This is clearly not the case 
for the 0.05 isochrones, with Worthey's (1994) being ~ 
0.5 magnitude redder in V-K. For metal-rich populations 



Table 1. V-K colours for simple stellar populations of age tQ = 
12 Gyr and metallicity Z. BBCFN94=Bertelli et al. (1994)-|-low 
mass extensions; W94=Worthey (1994); •=Worthey (1995). 
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such as giant ellipticals, it is imperative to be aware of, and 
understand the origin of, this discrepancy. 

Overlaying the Z=0.02 and Z=0.05 isochrones (only up 
to the red giant branch tip, for clarity), we see the difference 
quite graphically in Figure 2. The main sequences, sub-giant 
branches, and lower red giant branches agree for Z=0.05, but 
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V-K 

Figure 2. Solar and supersolar metallicity isochrones of age 12 
Gyr. The post-red giant branch components have been eliminated 
for clarity. W94=Worthey 1994; BBCFN94=BerteUi et al. 1994. 

the upper giant branches, as seen in this Mboi-(V-K) plane, 
increasingly diverge. This divergence is not evident in the 
Z=0.02 isochrones, shown for comparison in the same plot. 

The origin of this 0.5 magnitude discrepancy for 
0.05 seems to lie, at least in part, with the colour (V-K) - 
effective temperature relation used by the different groups. 
Whereas Worthey (1994) uses a metallicity-dependent scale 
for Z>Zq fluxes, Bertelli et al. (1994) do not (see their Sec- 
tions 4(ii)-4(iv)). In Figure 3 we show the main sequence and 
red giant branch (V-K)-log Teff relations used for the solar 
(solid curves) and Z=0.05 (open and filled circles) isochrones 
from the two sources. Worthey's (1994) Z=0.05 giant branch 
can be seen to be ~ 0.3 0.6 magnitudes redder in V-K, 
for a given temperature, when compared with his Z=0.02 
giant branch. In contrast, the Bertelli et al. (1994) Z=0.05 
giant branch relation overlays precisely their Z=0.02 rela- 
tion, a partial contributor to their significantly bluer giant 
branch, when compared with Worthey (1994), for super- 
solar metallicities. We should stress that the gross offset 
in the RGB (V-K)-Teff plane for Z<Zq, between Worthey 
(1994) and Bertelli et al. (1994), is masked/compensated by 
the two groups' different post-main sequence lifetimes and 
GB tips. This "conspiracy" means that the Z<Zq SSP in- 
tegrated colours are always within ~ 0.1 mag (in V-K) of 
each other. Assuming this conspiracy holds for the Z=0.05 
isochrones, which admittedly may not be entirely correct, 
implies that if Bertelli et al. (1994) scaled their GB redward, 
by what appears to be the more appropriate ~ 0.3 0.6 
mag (V-K), then the "usual" < 0.1 mag difference in the 
groups' colours would be recovered. It appears that Bressan, 
Chiosi & Fagotto (1994), who themselves use the Bertelli 
et al. (1994) isochrones, seem to have modified the latter's 
Z=0.05 (V-K)-log Teff relation, as their SSP integrated V-K 
is ~ 4.3, as opposed to the ~ 3.7 found in the digital version 
of the isochrones (at ta = 12 Gyr). This is ~ 0.15 mag- 



nitudes redder than Worthey's (1994) equivalent isochrone 
(Table 1, and Worthey's (1994) Figure 35), but does appear 
to bring their late-time integrated Z=0.05 SSP colours more 
in line with those of Worthey (1994). This is not to imply 
that Worthey's (1994) approach is the "last word" on the 
subject; indeed, we are of the opinion that the stellar evolu- 
tion models adopted by Bertelli et al. (1994) (and the subse- 
quent Padua Group compilations) are superior to those used 
by Worthey (1994) (in particular, the reduced post-main se- 
quence lifetimes in the Padua compilations seem preferable 
- Chariot, Worthey & Bressan 1995; Worthey 1995). Con- 
versely, though, we believe that Worthey's (1994) spectral 
library (specifically, at low temperature, where his detailed 
colour-temperature analysis is most commendable), is pre- 
ferred over that used in the Bertelli et al. (1994) isochrones. 

The situation becomes further complicated when we ex- 
amine the updated Bertelli et al. (1994) and Bressan, Chiosi 
& Fagotto (1994) isochrones (see Table 2 of Tantalo et al. 
1995). Tantalo et al. (1995) have updated their (V-K)-log 
Teff relation, such that it now resembles very closely that 
used by Worthey (1994) (recall Figure 3). Intuitively, follow- 
ing the argument of the previous paragraph, this should lead 
to (V-K) colours which are ~ 0.5 redder than the Bertelli 
et al. (1994) isochrones, at late-time. (Aside: We specu- 
late that this is indeed what was done by Bressan, Chiosi 
& Fagotto (1994), as it meshes with the colours quoted in 
their Table 3). Surprisingly, though, the (V-K) colours in 
Tantalo et al. (1995) are actually ~ 0.1 mag bluer than 
Bertelli et al. (1994) for Z=0.02, and roughly unchanged for 
Z=0.05. Again, we can only speculate, but we feel that the 
possible reason for this is that Tantalo et al. (1995) adopt 
an increased post-main sequence mass loss parametrisation 
{rj = 0.45, as opposed to the conventional rj = 0.35 used by 
Bressan, Chiosi & Fagotto 1994), which could contribute to 
their bluer SSPs through the consequent reduction in the 
RBG and AGB lifetimes. 

Disentangling all the differences between the Bertelli 
et al. (1994), Bressan, Chiosi & Fagotto (1994), Tantalo 
et al. (1995), and Worthey (1994) isochrones is a difficult 
task. A very important first step in this process can be found 
in the excellent study by Chariot, Worthey & Bressan (1995) 
into some of the uncertainties in modeling stellar popula- 
tions, to which the reader is directed for more details. We 
simply wish to illustrate what appears to be a key differ- 
ence between the Z>Zq isochrones that was encountered, 
independently, during our analysis, in order to understand 
the differences which appear during the calculation of the 
present-day photo-chemical properties of our model ellipti- 
cals. 

In what follows, we will provide parallel photometric 
evolution analyses based upon both the Bertelli et al. (1994) 
and Worthey (1994) isochrones. As to which set is "prefer- 
able" is left to the reader's discretion. 

As a final note to implementing the Bertelli et al. (1994) 
isochrones, there is an error in the running tabulation of the 
infrared colours in the digital version of their tables, which 
will be corrected in future releases (Bertelli 1994). 
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Figure 3. Colour-effective temperature relationsliip adopted for 
tlie solar and supersolar metallicity isochrones of Figure 2. Only 
the main sequence (MS) and red giant branch (RGB) are shown. 
Other symbols as in Figure 2. 



2.2 Photometric evolution 

Consider, first, the flux F through the bandpass AA for a 
composite stellar population (CSP) of metallicity Zv (given 
by equation 2), at time ta, as the sum of a series of i simple 
stellar populations (SSPs), each of age ta—t, and metallicity 
Z{t,). For each SSP isochrone, we sum over the wg-s initial 
mass bins, each j"^ bin having a luminosity Lj, resulting in 
an integrated flux 
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The subscript "G-i" indicates the nearest isochrone of age 
ta — ti to the metallicity Z{t,). (I>{m) oc is the initial 

mass function, by mass, where a slope x = 1.35 represents 
the canonical Salpeter (1955) value, ipi't) is the star forma- 
tion rate. 

Worthey (1994) (Z < 0.05) and Worthey (1995) (Z > 
0.05) isochrones were kindly provided by Guy Worthey with 
a temporal resolution of 0.1 Gyr. ~0.5 Gyr is the typical 
spacing for the Bertelli et al. (1994) tables. This was found 
to be more than adequate for our purposes. We are con- 
cerned predominantly with isochrones of age 10 ^ 12 Gyr 
for this work, as we will be comparing photo-chemical pre- 
dictions against those seen in the present-day local elliptical 
population. The earliest phases of spectral evolution will not 
covered in this work. 

Following Arimoto & Yoshii (1986), we assume that 
stars born in the discrete interval ^-^/Zk-i Zk \JZhZh-\-i\ 
are tied to the isochrone of metallicity Zh- For the 
Worthey (1994) isochrones, the subscript k runs from one to 
eleven, with Zi=0.0002, Z2=0.0005, Z3=0.0017, Z4=0.0053, 



Z5=0.0095, Z6=0.0169, Z7=0.0301, Z8=0.0361, Z9=0.0434, 
Zio =0.0534, and Zn =0.0950. The first and last intervals 
are Z=0 0.00030 and Z=0. 07127 oo, respectively. For 
the Bertelli et al. (1994) isochrones, the subscript i runs 
from one to six, with Zi=0.0004, Z2=0.0010, Z3=0.0040, 
Z4=0.0080, Z5=0.0200, and Z6=0.0500, with first and last 
intervals of Z=0 0.00063 and Z=0. 03162 oo, respec- 
tively. 

We calculate a V-band luminosity-weighted metallic- 
ity which in spirit is similar to Arimoto & Yoshii (1987), 
but differs in that we calculate the mean metallicity <Z>, 
whereas in their original work they determined < [Z] > 
(i.e. they replaced Z with [Z] in equation 2). On the sur- 
face this would seem to be a trivial difference, but as 
it turns out, this has crucial implications for estimating 
the metallicity for the model galaxy. Using the strict Ari- 
moto & Yoshii (1987) formalism results in any low metal- 
licity component being weighted far too heavily. A simple 
"back-of-the-envelope" example which illustrates this point 
might be useful - let us take a simple scenario whereby we 
have ten stars of equal luminosity, one star of which has 
Z=0.02 (i.e. [Z]=-|-0.000) and the other nine have Z=0. 00002 
(i.e. [Z]=-3.000). Working in "log" space, as Arimoto & 
Yoshii (1987) did, would yield a mean luminosity-weighted 
metallicity of < [Z] > =((-3.0x9)-|-(0.0xl))/10.0=-2.70, 
whereas strict adherence to equation 2 would lead to 
<Z> = ((0.00002x9)-|-(0.02xl))/10.0= 0.002018, or [<Z> 
]=-0.996. In other words, equation 2 would imply a mean 
metallicity 50 times higher than Arimoto & Yoshii (1987) 
would find. 

It is clear from the work of Worthey (1994,1995) that 
using equation 2 is more representative of "reality" as most 
spectral indices scale much closer to Z (i.e. number of ab- 
sorbers) than log Z. Of course, modelling the actual absorp- 
tion features (e.g. Mg2 line index) is the full and proper 
method for assigning metallicities, as the dependence upon 
Z is neither perfectly linear nor logarithmic, but for our pur- 
poses working with the linear regime for our "photometric" 
metallicity determinations is more than adequate. 
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It is important to be aware that it is common prac- 
tice to only present a mass-weighted metallicity (e.g. Mat- 
teucci 1994; Bressan, Chiosi & Fagotto 1994; Elbaz, Arnaud 
& Vangioni-Flam 1995), which may overestimate the actual 
observed metallicity by > 0.2 dex, even for massive galax- 
ies (e.g. Yoshii & Arimoto 1987; Gibson 1995). Comparison 
with the observed metallicity-luminosity relation should be 
made using the luminosity- weighted (i.e. equation 2), rather 
than the mass-weighted, metallicity. 



2.3 The bimodal star formation model 

The general bimodal star formation model for the evolu- 
tion of elliptical galaxies was introduced by Arnaud et al. 
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(1992), and described more quantitatively by Elbaz, Arnaud 
& Vangioni-Flam (1995). As this particular formalism for 
"bimodal" evolution is a popular one (e.g. its use in mod- 
eling the ultraluminous IRAS 10214+4724 - Elbaz et al. 
1992), we shall use their nomenclature, and refer the reader 
to their papers for any details not apparent in the descrip- 
tion which follows. 

The evolution of an elliptical is considered in four dis- 
tinct phases. Phases 1 and 111 are the two in which star 
formation occurs. The first phase is strictly a high mass 
"burst" mode in which a Salpeter (1955) IMF is assumed, 
but with a high low-mass cut-off of mi = 3.0 Mq (and nor- 
mal upper mass of 100 Mq), similar to that seen in certain 
starburst phenomena (Chariot et al. 1993). The star for- 
mation is intense, and short-lived (~ 15 25 Myrs for the 
initial masses considered below), due to the ensuing large 
SNe Type 11 rate, which subsequently drives a galactic wind 
for ~ 30 ^ 45 Myrs (i.e. the so-called Phase 11), thereby 
halting temporarily further star formation. Finally, the sys- 
tem enters a longer-lived (~ 1.5 2.5 Gyrs) phase (i.e. 
Phase 111) of milder star formation, again with a Salpeter 
(1955), but this time over the full range of masses observed 
today (i.e. ~ 0.1 100 Mq). Phase 111 star formation ter- 
minates with the onset of a second, smaller, galactic wind, 
once the thermal energy in the ISM exceeds that of its bind- 
ing energy (Larson 1974). The subsequent, quiescent. Phase 
IV lasts until the current epoch. 

Star formation during Phase 1 is exponential in time (t 
in Gyr): 



^i(t) = j.ie-'/^Mg(0) [MQ/Gyr] 



(3) 



where Mg{0) is the initial gas mass in Mq; the astration 
parameter in Phase 1 is taken to be = 40 Gyr~^ ; and, ipi 
decays exponentially with r = 0.05 Gyr. The gas mass and 
mass of metals Z in the gas during Phase 1 evolve as 



dt 



-Mi)+ J (t>im)ip{t-Trr^)R{m)dm (4) 

max[mf ,3.0] 



and 

dMzjt) 
dt 



(j>{m) 



ip(t-Tm)m1l ^dm(5) 



max[mf ,3.0] 



respectively. The second terms on the right-hand sides of 
equations 4 and 5 represent the ejection rate of gas and 
metals from dying stars. The Schaller et al. (1992) Z=0.02 
stellar lifetimes Tm were adopted. Parallel models run with 
the Schaller et al. (1992) Z=0.001 lifetimes were negligibly 
different, corroborating the conclusion of Gibson (1995) that 
the lack of full stellar lifetime metallicity-dependence is not 
a major source of uncertainty. The yields „ come from 
the recent metallicity-dependent compilation of Woosley & 
Weaver (1995), and we use the "mass in/mass out" formal- 
ism of Timmes, Woosley & Weaver (1995) (as opposed to 
Talbot & Arnett's (1973) matrix form) for solving equation 
5. R{m) = 1 — Wm/m represents the fractional mass of a 
star of initial mass m and remnant mass Wm (from Prant- 
zos, Casse & Vangioni-Flam 1993) which is ejected back into 
the ISM after its lifetime Tm- i/'("^) is in the IMF, 

by mass, and x = 1.35 is the classical Salpeter (1955) slope 



used here. We stress that the ultimate conclusions of Section 
3, regarding the bimodal star formation model, are not de- 
pendent upon the adopted stellar lifetimes, remnant masses, 
or stellar yields. 

In Phase 11, star formation is taken to cease for the 
duration of the galactic wind, and the total mass in the 
system Mg is assumed to evolve (i.e. decrease due to the 
steady wind outflow) as 



-aMsit) [MQ/Gyr] 



(6) 



dt 

where Mg = Mg{0) — Mgj, and Mgj is the total mass ejected 
up to time t; is the characteristic exponential timescale 
for mass loss - a = 18 Gyr— 1. The mass of metals in the 
remaining gas at each step is also adjusted accordingly. 

The outflow stops once the gas has cooled to a level 
below that of the gas' binding energy, and the Phase 111 
star formation regime is entered, with 



fiu{t) = ,,2M,{t) [MQ/Gyr] 



(7) 



In other words, star formation is taken to be proportional to 
the available gas mass, but at a significantly reduced level to 
that encountered in Phase 1. The astration parameter 1^2 = 2 
Gyr-i. 

The evolution of the gas and metal mass in Phase 111 
parallels that of equations 4 and 5, except that ipi is obvi- 
ously replaced with ^iii, and the lower limits on the integrals 
are replaced with max[mt, 0.1]. During Phase 111 there is also 
the additional complication of gas being returned to the sys- 
tem from stars which were born during Phase 1, but whose 
ejecta did not partake of the outflow of Phase 11. Thus for 
the first ~ 0.3 Gyr (i.e. the lifetime of a 3 Mq star minus the 
duration of the outflow) of Phase 111 there is an additional 
"ejecta" term for the dying lower-mass Phase 1 stars. 

To ensure compatibility with the Elbaz, Arnaud 
& Vangioni-Flam (1995) results, we have set explicitly 
the Phase 1, 11, and 111 demarcation (e.g. wind initia- 
tion/completion) times given in their Table 2b, as opposed 
to recalculating explicitly the gaseous thermal energy evolu- 
tion. While the numerical examples highlighted in this paper 
were selected on the basis of their being the favoured Elbaz, 
Arnaud & Vangioni-Flam (1995) models, the models hold, 
in general, for most permutations of this simple bimodal sce- 
nario. Let us now turn to the specific photo-chemical pre- 
dictions for the model elliptical galaxies. 



3 DISCUSSION 

We ran models identical to Elbaz, Arnaud & Vangioni- 
Flam's (1995) preferred grid (Models A^D in their Table 
2a). The initial gas masses were 5 X 10^ 5 X 10^", 5 X 10", 
and 5 X 10 Mq. No dark matter halos were considered for 
these particular models, although their Models E^H, with 
said halos, show this to be an unimportant consideration for 
what follows. 

Figure 4 shows the evolution of each galaxy's gas mass 
as a function of time. The gas mass decreases during Phase 1 
due to star formation, and continues to decrease in an over- 
all sense during Phase 11 (due to the SNe Type 11-driven 
outflow). The vertical line marks the onset of Phase 111. 
We draw attention to the onset of Phase 111 because only 
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Figure 4. Evolution of the gas mass fraction of the ISM for the 
bimodal star formation models of Section 2.3. Only stars formed 
during Phase III are still alive at the current epoch. 

stars formed during this phase contribute directly to the ob- 
served photo-chemical properties at the current epoch. The 
gas mass increases slightly for a few tenths of a Gigayear 
after the initiation of Phase III due to to gas returned by 
dying stars formed in Phase I. Once the last of these (m = 3 
Mq) have died, the gas mass decreases drastically once again 
due to the star formation governed by equation 7. The be- 
haviour of the gas mass fraction evolution is very similar for 
each of the models concerned. Figure 4 is directly compa- 
rable with the ISM curve in Figure 3 of Elbaz, Arnaud & 
Vangioni-Flam (1995), although it should be noted that we 
have normalised to the initial gas mass whereas they nor- 
malise to the current galactic mass. The single curve that 
they show (Mg(0) = 5 x 10"'"' Mq) is indistinguishable from 
that shown in our Figure 4, once this difference in normali- 
sation is accounted for. 

Of more interest are the results shown in Figure 5. Here 
we see the evolution of the global metallicity of the ISM for 
the four different initial masses. The rapid rise of Z during 
Phase I is traced to the high initial star formation rates, and 
more importantly, the imposed high-mass mode of the IMF. 
The ISM metallicity passes the solar value in ~ 14 Myrs. 

The metallicity during the initial few tenths of a Gi- 
gayear in Phase III increases very slowly, or even decreases, 
from its initial value, due to dilution effects from lower mass 
stars formed at lower Z in Phase I dying during the early 
portion of Phase III. This effect is enough to compensate 
for the increased metal production from the high-mass stars 
now dying in Phase III which themselves were formed in 
the pre-enriched (from Phase I stars) Phase III ISM. After 
t Si 0.3 Gyr, the global Z increases dramatically for ~ 1 Gyr 
due to the cessation of this Phase I dilution influence. 

Figure 6 shows the metallicity evolution as measured 
by [Z] and [Fe], for the Mg{0) = 5 x lO" Mq model. Allow- 
ing for the different normalisation, the [Fe] curve is indis- 
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Figure 5. Evolution of the ISM metallicity Zg for the bimodal 
star formation models of Section 2.3. Only stars formed during 
Phase 111 are still alive at the current epoch. 

tinguishable from that shown in Figure 4 of Elbaz, Arnaud 
& Vangioni-Flam (1995) (note: they normalise to the total 
initial mass). We note, further, that [Fe]si[Fe/II] for these 
models. [Fe/H] roughly tracks the global metallicity [Z] in 
Phase III, although offset by typically ~ 0.2 0.3 dex (al- 
though always super-solar in Phase III, regardless of metal- 
licity "definition"). See also the comment at the end of this 
section regarding the SNe la formalism adopted by Elbaz, 
Arnaud & Vangioni-Flam (1995) which may mean that the 
iron abundance of Figure 6 is actually an underestimate. 

In this work we generally refer to the evolution of the 
global metallicity Z, as opposed to the metallicity denoted 
by [Fe/H], simply because the observed "[Fe/H]" tabulated 
in Terlevich et al. (1981), for example, are actually closer 
to [Z] than to [Fe/H], per se (Terlevich 1995), as they are 
indirectly assigned from a combination of spectroscopic indi- 
cators (Mg2, in this case) and photometric calibrators (e.g. 
Barbuy 1994, and references therein). This is confirmed by 
Gonzalez & Gorgas (1995). Obviously, neither [Z] nor [Fe/H] 
are ideal, but as the current version of our code is only "pho- 
tometric" in nature (i.e. no spectral line indices), we chose 
the best compromise and compared our V-band luminosity- 
weighted metallicity [<Z>]v against the observations shown 
in Figure 7. Elbaz, Arnaud & Vangioni-Flam (1995) only 
report the mass- weighted [Fe/H] metallicity for their mod- 
els, so this difference should be kept in mind, although we 
stress that any differences which might crop up due to this 
[Z] versus [Fe/H] metallicity "definition" do not affect the 
conclusions of this paper. 

Recalling that only stars formed during the "normal" 
IMF mode of Phase III are still alive today to contribute to 
the elliptical's luminosity (and hence photo-chemical prop- 
erties), we can anticipate immediately from Figure 5 that 
there is going to be a problem with their bimodal star for- 
mation scenario. Specifically, each of the Phase III curves 
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Figure 6. Evolution of the global ISM metallicity [Z] and its iron 
content [Fe] for the bimodal star formation model of initial mass 
5 X 10"^"^ Mq. The [Fe] curve matches that of Elbaz, Arnaud &l 
Vangioni- Flam's (1995) Figure 4. 

only span < 0.2 dex in [Z], whereas observations show a 
metallicity spread of ~ 3.5 dex (e.g. Minniti et al. 1995, 
and references therein). Not only is their implied metallic- 
ity distribution virtually non-existent, but there are no stars 
with Z< (2 3)Zq (again, at odds with the inferred low- 
metallicity components in ellipticals - e.g. Rose 1985). We 
conclude that just on the basis of this rough "chemical ar- 
gument" there would seem to be concern for the bimodal 
star formation scenario. Still, let us play devil's advocate 
and claim that since there is only minimal direct evidence 
for metallicity spreads in ellipticals (e.g. a dispersion > 2 
dex is seen in the giants in the nearby dwarf ellipticals M32 
and NGC 147 - Davidge & Jones 1992; Davidge 1994 - the 
bulk of the evidence is indirect in nature - i.e. based upon 
spiral bulge stars - e.g. Minniti et al. 1995), this chemi- 
cal argument alone is not enough for dismissing the simple 
bimodal model. Unfortunately, the problems become even 
more acute when we look at their predicted present-day 
colour-metallicity-luminosity relations - i.e. direct observa- 
tional tests which do hold for local ellipticals. 

The framework adopted by Elbaz, Arnaud & Vangioni- 
Flam (1995) (and most of the wind model papers referenced 
thus far) - coeval galaxies with the present-day metallicity- 
luminosity relation driven by the cessation of star formation 
due to the galactic winds - is a typical, and elegantly sim- 
ple, one in which to work, but it should be stressed that this 
may not be a unique interpretation. For example, Worthey, 
Trager & Faber (1995) speculate that the observed rela- 
tions actually reflect a luminosity-weighted age-metallicity 
relation with the most metal-rich ellipticals being ~ 2 Gyr 
old (i.e. later phases of star formation taking place in pre- 
enriched gas from an older epoch of star formation), and the 
most metal-poor ~ 15 Gyr (i.e. little or no late-time star 
formation). These Worthey, Trager & Faber (1995) ages are 



based on C24668 line strengths and an average of Balmer 
indices, but as they themselves stress, it only takes a "frost- 
ing" of a fraction of a percent of a very young population 
on top of a basically old population to impact dramatically 
upon the Balmer line strengths, making a galaxy appear 
much younger than a mass-weighted mean age (Section 3 
of Worthey, Trager & Faber 1995). Obviously, adopting any 
sort of "mean" age or "mean" metallicity will mask, at some 
level, the true complexity of the star formation process in 
ellipticals, a model limitation which should always be kept 
in mind. 

Figure 7 shows the predicted metallicity-luminosity re- 
lation for the grid of models described above (the two dotted 
lines represent the relevant photo-chemical predictions de- 
rived from the two isochrone sources - in this plane, the 
Worthey (1994) and Bertelli et al. (1994) predictions are 
virtually indistinguishable). The solid curve represents the 
template of models outlined in Gibson & Matteucci (1995) . 
As we could anticipate from Figure 5, the Elbaz, Arnaud 
& Vangioni-Flam (1995) present-day luminosity-weighted 

* The Gibson &l Matteucci (1995) template adopts a single, some- 
what flatter-than-Salpeter (1955) IMF {x « 1), with mi = 0.2 
Mq and mjj = 65 Mq; Worthey 's (1994) isochrones were as- 
sumed. Early galactic wind times (i.e. ^gw ^ 0-1 Gyr), similar 
to those favoured by Bressan, Chiosi &z Fagotto (1994) (whose 
colour-luminosity predictions are similar to the Gibson &z Mat- 
teucci (1995) grid's in Figure 8), were found, despite our con- 
servative approach to pre-SN stellar wind energy deposition to 
the ISM (Gibson 1994a), and their extreme view that SNe are 
only negligible contributors in setting tow • Indeed, despite the 
analysis of Gibson (1994a), Bressan, Chiosi &z Tantalo (1995) 
and Tantalo et al. (1995) remain adamant that if the colour- 
metallicity relation of ellipticals is to be understood, stellar wind 
energy must dominate. How then can the Gibson &z Matteucci 
(1995) photo-chemical predictions (and, in particular, tow) be 
comparable with those of Bressan, Chiosi &z Fagotto's (1994)? The 
answer lies in the treatment of SN remnant dynamics and ther- 
mal energetics. Bressan, Chiosi &z Fagotto's (1994) treatment (and 
indeed, that of Arimoto &z Yoshii 1987; Matteucci &z Tornambe 
1987; Angeletti &z Giannone 1990; to name but a few) is the clas- 
sic "evolution in isolation" formalism in which remnants are as- 
sumed to expand and cool in complete isolation ad infinitum, 
resulting in the inefficient transfer of thermal energy to the ISM. 
This fundamental assumption inherent in most previous galactic 
wind models is incorrect, and does not reflect the complexity of 
reality — in reality, as recognised by Larson (1974), Dekel &z Silk 
(1986), Babul & Rees (1992), Gibson (1994b), and Nath & Chiba 
(1995), shells quickly come into contact and overlap with neigh- 
bouring, expanding shells, affecting strongly both their dynamics 
and energetics. Subsequent SNe explosions occur in the resultant, 
expanding, rarefied bubble — i.e. "superbubble" (Tomisaka 1992) 
— with very little adiabatic or radiative cooling energy loss. In 
reality then, the SNe energy transfer to the ISM is very efficient, 
and early galactic wind times are a natural outcome; there is no 
need to resort to what we feel are physically implausible assump- 
tions regarding the stellar wind energy, as was done by Bressan, 
Chiosi & Fagotto (1994), Bressan, Chiosi & Tantalo (1995), and 
Tantalo et al. (1995), in order to fit the galactic colour-luminosity 
relation. The exact formalism adopted for the SNe energetics is 
given by Model B^ of Gibson (1994b), whereas that adopted by 
Bressan, Chiosi & Fagotto (1994) would correspond to Gibson's ( 
1994b) Model Aq . Further details can be found in Gibson (1994b) , 
Gibson & Matteucci (1995), and Gibson (1995). The latter work 
addresses the Bressan, Chiosi &z Fagotto (1994) paper specifically. 
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Figure 7. Observed metallicity-lmninosity relation for dwarf 
(open circles: Smith 1985) and "normal" (open squares: Sil'chenko 
1994 — filled circles: Terlevich et al. 1981) ellipticals. The 
[Mg/Fe] =+0.25, Mg2-[Z] calibration of Barbuy (1994) has been 
adopted for the normal/giant ellipticals. The dotted lines repre- 
sent the predictions for the bimodal star formation/IMF models of 
Elbaz, Arnaud &z Vangioni-Flam (1995) (the two predictions come 
from the Worthey (1994) and Bertelli et al. (1994) isochrones, 
and are virtually indistinguishable), while the solid curve shows 
the template models of Gibson &z Matteucci (1995). 
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Figure 8. V-K colour- luminosity relation for Virgo (open circles) 
and Coma (filled circles) cluster ellipticals and lenticulars (Bower, 
Lucey & EUis 1992). A Virgo distance modulus of {V - My)o = 
31.54 is assumed, and a shift of A[V — My)o = 3.58 has been 
applied to the Coma sample. Data for the dwarfs from Thuan 
(1985). The dotted lines represent the predictions for the bimodal 
star formation/IMF models of Elbaz, Arnaud &z Vangioni-Flam 
(1995) (upper line implied use of the Worthey (1994) isochrones, 
while the lower line used the Bertelli et al. (1994) isochrones), and 
the solid curve shows the template models of Gibson &z Matteucci 
(1995). 



metallicities are considerably greater than what is seen in 
nearby ellipticals - < [Z] >v is discrepant with the observa- 
tions by ~ 0.2 1.0 dex, the predictions getting progres- 
sively worse for lower masses. Again, as anticipated from 
Figure 7, it is not difficult to trace the problem to the com- 
plete lack of low-metallicity stars at ta = 12 Gyr. The high 
mass star formation mode of Phase I may be an efficient 
mechanism for ejecting vast quantities of O and Fe to the 
cluster ICM, but its unavoidable enrichment of the Phase III 
ISM means that no sub-solar metallicity stars remain today. 
This is at odds with the direct and indirect arguments for 
sub-solar metallicity components in nearby ellipticals, put 
forth by Rose (1985), Arimoto & Yoshii (1987), Hardy et al. 
(1994), and Minniti et al. (1995), amongst others. 

The problems are equally apparent in the colour- 
luminosity plane. In Figure 8 we show once again, the pre- 
dicted V-K for the Elbaz, Arnaud & Vangioni-Flam (1995) 
bimodal star formation models (dotted lines, where the 
upper line shows the colours predicted from the Worthey 
(1994) isochrones, and the lower line, the Bertelli et al. 
(1994) isochrones), as well as those for the proposed tem- 
plate of Gibson & Matteucci (1995) (solid curve). Again, 
the Gibson & Matteucci (1995) models use the more conven- 
tional Arimoto & Yoshii (1987) IMF, and satisfy the obser- 
vational constraints of both the ICM and the underlying stel- 
lar populations. The Elbaz, Arnaud & Vangioni-Flam (1995) 
models are far too red - by ~ 0.8 1.4 mag, when using 
the Worthey (1994) isochrones, and by ~ 0.4 1.0 mag. 



even when using the bluer Bertelli et al. (1994) isochrones. 
Recalling that there are no Z< 0.03 stars left at the cur- 
rent epoch in their models, this should not be too surpris- 
ing. Table 1 shows us that the bluest component in these 
CSPs would then only have V-Ksi3.9 (Worthey) or Si3.6 
(Bertelli et al. ). We remind the reader that the Gibson & 
Matteucci (1995) template was tuned to fit the observed the 
colour-metallicity-luminosity relations of Figures 7 and 8 via 
use of the Worthey (1994) isochrones. Adopting the Bertelli 
et al. (1994) isochrones, in their place, would necessitate the 
change of some of the specific details (e.g. star formation ef- 
ficiency J/, galactic wind epoch taw - see also Gibson 1995) 
in order to recover the observed relations, but the general 
conclusion of the current study would remain unchanged. 

We mention in passing that the preferred grid of mod- 
els from Elbaz, Arnaud & Vangioni-Flam (1995) (Models 
A^ D) assumed a Type la SN rate formalism which is un- 
evolving in time, given roughly by the present-day observed 
rates (Turatto, Cappellaro & Benetti 1994). Considering 
that adopting more physically significant models for SNe 
la evolution (e.g. Figure 1 of Matteucci & Tornambe 1987) 
implies that said rate in the initial ~ 2 Gyrs of star forma- 
tion should be ~ 20 times greater than the present-day value 
(for the singly-degenerate Type la SNe progenitor model of 
Whelan & Iben 1973), there is the potential to have under- 
estimated the ISM iron content seen in Phase III of Figure 6. 
This would have the tendency to drive down their resultant 
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stellar [Mg/Fe] predictions (which currently are consistent 
with the observations of Worthey, Faber & Gonzalez 1992), 
and increase their stellar [Fe/H], thereby worsening the al- 
ready poor colour-metallicity-luminosity predictions shown 
in Figures 7 and 8. 



4 SUMMARY 

We re-iterate the conclusions drawn in Section 3 - the 
simplest bimodal IMF/star formation models for elliptical 
galaxies, with galactic winds (i.e. those of Arnaud et al. 
(1992) and Elbaz, Arnaud & Vangioni-Flam (1995)), may 
satisfy the ICM observational constraints, but only at the ex- 
pense of the resultant predicted present-day photo-chemical 
properties for the underlying ellipticals. Specifically, the pre- 
dicted V-K colours are ~ 0.5 1.4 magnitudes too red, the 
luminosity-weighted metallicities ~ 0.2 1.0 dex too high, 
and the intrinsic stellar metallicity dispersion is virtually 
non-existent (i.e. < 0.2 dex, versus the > 3.5 dex implied 
observationally). Renzini et al. (1993) have also criticised 
the simple bimodal models, albeit on more general terms 
than those presented in this paper. We conclude that the 
more conventional models with flatter-than-Salpeter (1955) 
IMFs and efficient SN thermal energy transfer to the ICM, 
and which use the photometric techniques outlined in this 
paper, are still the preferred scenario for elliptical galaxies 
with galactic winds (Gibson & Matteucci 1995), in our opin- 
ion. 
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